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The performance of nano-scale energy conversion devices is studied in the framework of state- 
space models where a device is described by a graph comprising states and transitions between them 
represented by nodes and links, respectively. Particular segments of this network represent input 
(driving) and output processes whose properly chosen flux ratio provides the energy conversion 
efficiency. Simple cyclical graphs yield Carnot efficiency for the maximum conversion yield. We give 
general proof that opening a link that separate between the two driving segments always leads to 
reduced efficiency. We illustrate this general result with a simple model of an organic photovoltaic 
cell, where such an intersecting link corresponds to non-radiative carriers recombination and where 
the reduced maximum efhciency is manifested as a smaller open-circuit voltage. 

PACS numbers: 05.70.Ln, 05.60.-k,73.50.Pz 


Modeling of molecular scale engines such as photo¬ 
voltaic cells, thermoelectric devices and light emitting 
diods is often done in terms of a state space description 
[l[. In such state-space models the processes underlying 
the device operation are described as transitions between 
microscopic system states, leading to the engine descrip¬ 
tion in terms of a graph comprising nodes (=states) con¬ 
nected by bonds that represent transitions between them. 
On this network, the time evolution is determined by the 
rates associated with each bonds and fluxes associated 
with the non-equilibrium dynamics flow along intercon¬ 
nected linear and cyclical paths. Particular segments of 
this network represent the input (driving) and output 
processes, where the corresponding rates deviate from re¬ 
strictions imposed by equilibrium thermodynamics, thus 
maintaining the system in a non-equilibrium state. Sim¬ 
ilar models are often used to describe and analyze bio¬ 
chemical networks 

An example is shown in Fig. [T] that depicts a single cy¬ 
cle model (henceforth referred to as model SC) for a pho¬ 
tovoltaic cell comprising N states that represent a charge 
carrier occupying different sites in the system and one 
state (denoted 0) in which the charge carrier is absent. 
Transition rates between system states j = 0,1, • • • ,N, 
are determined by state energies Ej, the bias voltage 
V and the local temperatures. The latter are taken 
equal to ambient temperature T everywhere except at 
the s ^ s J- 1 transition that is driven by the sun tem¬ 
perature Ts- The different rates satisfy detailed balance 
conditions: 





(b) - 



s s+1 



N 


FIG. 1. (color online), (a) A simple graph (model SC) rep¬ 
resenting an ideal (no non-radiative losses) photovoltaic cell. 
The system is driven away from equilibrium by the processes 
that take place at the red and blue segments, (b) Dynamics 
at the red segment: States s and s -f 1 are ground and ex¬ 
cited states of a dye molecule and transition between them 
is driven by the “sun temperature” Ts- (c) Dynamics at the 
blue segment. The transitions 0^1 and 0 ^ A represent 
changes in the number of electrons on the molecule due to its 
coupling to the left and right electrodes, respectively. Non¬ 
equilibrium is imposed by the voltage difference V between 
the two electrodes, prj — pi = E (the electron charge e is set 
to 1 here and below). 
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where kij is the rate coefficient for the j ^ i transition 
and /3 = {ksT^^ , /3s = {kBTs)~^- The cycle affinity A 
is defined by [5| 


A = -\n]C 


( 4 ) 


where /C is the ratio between products of forward and 
backward rates 
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The open-circuit (OC) voltage is determined by the 
condition that the drivings associated with the voltage 
V and with the difference between T and Ts balance 
each other so that the steady state current through the 
wstem vanishes. In this case, the cycle affinity vanishes 
[6|, namely K, — 1. Denoting the corresponding voltage 
^J'R - ^J'L by Voc and using Eqs. (HD-® this leads to 

= A£ = £.«-£;.. (6) 

The l.h.s is limj_s.o ??(J) where r]{J) = V{J)/AE is the 
thermodynamic (“internal”) efficiency of the device and 
J is the steady state current (same through all links). 
Equation m thus yields the Carnot expression rjc for 
the efficiency in the zero power limit. It also highlights 
the OC voltage as an important attribute in the quest 
to understand cell performance In Ref. [Tl| we 

have derived two important generalizations of this result: 
First, if the transition s ^ s-l-1 results from the combina¬ 
tions of radiative and non-radiative processes whose for¬ 
ward/backward ratios are determined by temperatures 
Ts and T, respectively, Eq. ([6]) is replaced by a simi¬ 
lar expression in which Ts is substituted by an effective 
temperature [l^ 
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where p = is the ratio between the non- 

radiative and radiative relaxation rates. Second, if one of 
the cycle processes involves charge separation that comes 
at an energy cost (exciton binding energy) Ts, Eq. ([6]) 
is replaced by 0 
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FIG. 2. (color online). Models characterized by intersected 
graphs - an new pathway is added between sites m and n. 
(a) The intersecting pathway does not separate the two driv¬ 
ing (red and blue) segments, (b) The bond n — m separates 
between the two driving segments which now sit on different 
subcycles. 


The later situation characterizes the operation of organic 
photovoltaic cells 1,1-T^. 


Model SC and Eqs. ®-® constitute a convenient 
framework for analyzing the OC voltage of photovoltaic 
cells, and, as shown in Ref. can be used as a starting 
point for studying optimal performance in finite power 
operations, provided that information on the system elec¬ 
tronic structure and dynamics can be incorporated into 
such state-space model. However, this single cycle model 
is oversimplified; realistic state-space models contain in¬ 
tersecting pathways that makes the analysis that leads 
to Eqs. ®-® more involved. Indeed, in state-space 
models of photovoltaic cells such intersecting pathways 
represent non-radiative losses, e.g., electron-hole recom¬ 
bination processes that are expected to reduce efficiency. 

As graphical structures, such intersecting pathways 
divide the cycle of model SC into separate cycles (see 
Fig. [2D, making it possible to have steady states in which 
the current in the power-extracting segment vanishes 
while internal currents exist elsewhere in the system. 
While intuitively we expect lower efficiencies in such sit¬ 
uations, no general proof for this is available. 


In this Letter we consider two generalizations of model 
SC shown in Fig. |21 In one [Fig. HJa)], an intersect- 
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ing pathway divides the system cycle into two subcycles 
without separating between the two driving segments. 
In the second [Fig. [5Jb)], the driving segments are posi¬ 
tioned on different cycles. We show that (a) Intersecting 
pathways of the type shown in Fig.[21[a) do not affect the 
OC voltage, hence their zero power operation is charac¬ 
terized by the Carnot efficiency, (b) Intersecting path¬ 
ways of the type shown in Fig. [2l[b) lead to zero power 
efficiency smaller that the Carnot value, manifested for 
photovoltaic cells be a lower OC voltage. 

We start with Model SC and take note of the two sites, 
m and n, between which a new pathway will be enabled. 
The rate ratio, Eq. dSD, can obviously be written as K. = 
K- 1 IC 2 where /Ci = • • • kYn+l^m/kn—l^n * * * 

and /C 2 — ^? 7 i,m— 1 ' ■' ^n-t-i,n/^m—i,m ' ‘' are ratios 

of rate products that correspond to the arcs I and 2 in 
Fig.H Now enable the pathway between sites n and m, 
that is, make kn^m and km,n different from zero. The 
graph is now divided into two subcycles, I and 2 in addi¬ 
tion to the original cycle. The corresponding cycle rate- 
product ratios are 

^ km,n p- kn^m 

Ni — Ni --; 1^2 — 1^2 T - 

^m,n 

so that /C 1 /C 2 = iCiiC 2 - 

Consider next the situation shown in Fig. EJa), where 
the two driving segments are on the same subcycle 1. 
Starting from the model SC in the OC situation so 
that 1C = IC 1 IC 2 = 1, enabling the n — m link does 
not lead to current on this bond. To show this note 
that since all rates on arc 2 satisfy Eq. o, we have 
K ,2 = = km,n/kn,m- Equation dH) then 

yields A ^2 = 1 and therefore ICi = 1. Consequently, con¬ 
necting the sites n and m will not create any extra current 
and subsequently does not affect the open-circuit voltage. 

In the situation of Fig. EJb), each of the subcycles cre¬ 
ated by enabling the n — m pathway contains a non¬ 
equilibrium element: driving by the sun in subcycle I 
and the driving by the voltage in subcycle 2. Focusing 
on the latter and using Eqs. m and m yields 

jC2 = ^ ( 9 ) 

kn^m 

where A/i = — ^il- Equation (|S]) then yields IC 2 = 

exp[—/3A^]. This indicates that with the m — n link en¬ 
abled the system cannot establish a state where currents 
through all links in subcycle 2 vanish unless A/x = 0. 

To understand the consequence for the open-circuit 
voltage, let us set and km,n = '^km,n and 

evaluate the derivative of the OC voltage, Vbc) with re¬ 
spect to A at A = 0. Foe is the value of Afi = fJ-R — ^J-L 
for which the current between the electrodes, (the blue 
segment in Eigs. [T] and [2]), therefore across all links in 
arc I, vanishes; namely kjj+iPj+i = kj^ijPj for all this 
links. From Eqs. m and m it follows that under OC 



FIG. 3. (color online). A 4-level spinless model of an or¬ 
ganic photovoltaic cell [Til. The states are: 0 - no excess 
electron, 1 - excess electron in the higher occupied molecular 
orbital (HOMO) of the donor, 2 - Following an optical transi¬ 
tion, electron is on the otherwise lower unoccupied molecular 
orbital (LUMO) of the donor (an exciton state). 3: Follow¬ 
ing charge separation, the electron has moved to the acceptor 
LUMO. The transitions 3^0 and 0 1 include electron 

transfer to the right and left electrode, respectively, 1^2 
is the optical excitation process and 2 3 is the exciton 

dissociation process (see Ref. [TJ for further details). The pro¬ 
cess 3^1 that divides the cycle to two subcycles represents 
radiationless recombination. 


conditions, the probabilities to find the system in states 
n and m are related to Fqc by 


_^ _ „ — P(VoC+Em—En) 
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where c is a positive number and Pj = Pje~^^C The 
derivative on the r.h.s is negative if {Pm > Pn)\=o, 
i.e. Foe < 0, and negative in the opposite case, 
{Pm < Pn)\=o', IAqc > 0 (this Can be inferred from 
the current expression, Jni — m — A {kn^mPm km,7iPn'j — 
{Pm “ -Pn) ■ Hence 


/ ^|Vbc| 

^ dX 


< 0 . 

A=0 
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Enabling the m — n link in the scheme of Eig. EJb) thus 
results in reduction of the OC voltage, namelt of the 
maximum thermodynamic efficiency. 

To illustrate this general observation for a particular 
system, we consider a simple 4-level model for an or¬ 
ganic photovoltaic cell comprising a donor acceptor (D- 
A) complex placed between two electrodes (see Fig. [3] 
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and Supplementary Material IS)- The D and A species 
are characterized by their HOMO levels, Edhi ^Ahj 
LUMO levels, £dl, £Al that translate into three relevant 
energy differences /S.E = £ 0 ^ — £dh (the optical gap), 
Ae = £dl ~ ^Al (interfacial LUMO-LUMO gap), and 
EDA,g = eAl ~£Dh (effective band gap). Charge carriers 
can enter and leave the system only by transitions be¬ 
tween the donor HOMO level and the electrode L, and 
between the acceptor LUMO level and the electrode R. 
These electrodes are represented by free-electron reser¬ 
voirs at chemical potentials {K = L;R). An impor¬ 
tant feature of this model is that a thermal fluctuation 
is needed to overcome the exciton binding energy Eb of 
an interfacial charge transfer exciton to generate a free 
charge carrier. With some simplifying assumptions, the 
system can be described by a kinetic scheme involving 
four states. A detailed description of this model and the 


underlying kinetic scheme was given in Ref. [11 


An additional link between states 1 and 3 correspond 
to non-radiative recombination that returns a free elec¬ 
tron on the acceptor back to the donor ground state. 
In the absence of this process, the open-circuit volt¬ 
age Vbc(A = 0) is found to be given by Eq. ([7]). The 
maximum thermodynamic efficiency is correspondingly 
rjc — Eb/AE, highlighting the loss associated with the 
exciton binding energy. 

RepeatiM the cycle analysis when the rates kis, k^i 
are finite 17| yields the following expression for the open- 
circuit voltage 


^Aic(A) = AErjc — Eb 
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setting fci 3 = Afci 3 and k^i = Xk^i. and taking the 
derivative with respect to A we find, using also Eqs.(S5- 
S8) 

/ ^Vbc \ ^ /;^r ^l3(fc32 + fcl2) A _ g-/3(Voo)>=0^ , 

\ dX J x=o ki2k23 V / 

(14) 

where (Vbc)A=o is given by ([7]). This is an explicit state¬ 
ment of the inequality m- 

In conclusion, we have presented a state-space network 
representation scheme of molecular scale engines and ap¬ 
plied it to models of photovoltaic cells. A general ar¬ 
gument based on cycle analysis yields conditions under 


which the opening of intersecting pathways in an other¬ 
wise cyclical graph leads to a decrease in the maximum 
efficiency. In the application to the operation of a pho¬ 
tovoltaic cell such intersecting link often corresponds to 
a carrier recombination process, and its consequence is 
manifested in a reduced open-circuit voltage. It should 
be emphasized that similar analysis can be carried out 
for more complex models that include, for example, po- 
laron formation and hot exciton dissociation, as well as 
for other types of nano-scale engines, provided that they 
can be modeled by kinetic transitions in the system state- 
space. As described elsewhere such an approach can 
be used for analysis of optimal performance under finite 
power operation, although with the unavoidable loss of 
the generality associated with a thermodynamic descrip¬ 
tion. 


[1] U. Seifert, Rep. Prog. Phys. 75, 126001 (2012). 

[2] T. L. Hill, Journal of Theoretical Biology 10, 442 (1966). 

[3] T. L. Hill, Free Energy Transduction and Bioehemieal 
Cycle Kinetics, 2nd ed. (Dover, New York, 2009). 

[4] E. Gerritsma and P. Gaspard, Biophysical Reviews and 
Letters 05 , 163 (2010). 

[5] J. Schnakenberg, Rev. Mod. Phys. 48, 571 (1976). 

[6] M. Einax and A. Nitzan, J. Phys. Chem. G 118, 27226 
(2014). 

[7] K. Vandewal, K. Tvingstedt, A. Gadisa, O. Inganas, and 
J. V. Manca, Nat. Mater. 8, 904 (2009). 

[8] T. Kirchartz, K. Taretto, and U. Ran, J. Phys. Ghem. 
C 113 , 17958 (2009). 

[9] W. J. Potscavage, A. Sharma, and B. Kippelen, Acc. 
Chem. Res. 42 , 1758 (2009). 

[10] P. K. Nayak, J. Bisquert, and D. Cahen, Adv. Mater. 
23 , 2870 (2011). 

[11] M. Einax, M. Dierl, and A. Nitzan, J. Phys. Chem. C 
115 , 21396 (2011). 

[12] This result is obtained from Eq. (26) of Ref@ by using 
the detailed balance relationships between the different 
rates. 

[13] N. C. Giebink, G. P. Wiederrecht, M. R. Wasielewski, 
and S. R. Forrest, Phys. Rev. B 83 , 195326 (2011). 

[14] L. J. A. Koster, S. E. Shaheen, and J. C. Hummelen, 
Adv. Energy Mater. 2 , 1246 (2012). 

[15] M. Gruber, J. Wagner, K. Klein, U. Hermann, A. Opitz, 
M. Stutzmann, and W. Briitting, Adv. Energy Mater. 
2 , 1100 ( 2012 ). 

[16] M. Einax, M. Dierl, P. R. Schiff, and A. Nitzan, Euro- 
phys. Lett. 104 , 40002 (2013). 

[17] See Supplemental Material at http://[.] 










